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ABSTRACT 

The recently initiated Arecibo Legacy Fast ALFA (ALFALFA) survey aims to map 
~ 7000 deg^ of the high galactic latitude sky visible from Arecibo, providing a HI 
line spectral database covering the redshift range between -1600 km and 18,000 
km s~^ with ~ 5 km resolution. Exploiting Arecibo's large collecting area and 
small beam size, ALFALFA is specifically designed to probe the faint end of the HI 
mass function in the local universe and will provide a census of HI in the surveyed 
sky area to faint flux limits, making it especially useful in synergy with wide area 
surveys conducted at other wavelengths. ALFALFA will also provide the basis for 
studies of the dynamics of galaxies within the Local and nearby superclusters, will 
allow measurement of the HI diameter function, and enable a first wide-area blind 
search for local HI tidal features, HI absorbers at z < 0.06 and OH megamasers in the 
redshift range 0.16 < z < 0.25. Although completion of the survey will require some 
five years, public access to the ALFALFA data and data products will be provided in 
a timely manner, thus allowing its application for studies beyond those targeted by 
the ALFALFA collaboration. ALFALFA adopts a two-pass, minimum intrusion, drift 
scan observing technique which samples the same region of sky at two separate epochs 
to aid in the discrimination of cosmic signals from noise and terrestrial interference. 
Survey simulations, which take into account large scale structure in the mass distribution 
and incorporate experience with the ALFA system gained from tests conducted during 
its commissioning phase, suggest that ALFALFA will detect on the order of 20,000 
extragalactic HI line sources out to z ~ 0.06, including several hundred with HI masses 
Mhi < 10^-5 Mq. 

Subject headings: galaxies: spiral; — galaxies: distances and redshifts — galaxies: halos 
— galaxies: luminosity function, mass function — galaxies: photometry — radio lines: 
galaxies 



1. Introduction 

The first 21 cm line detection of an extragalactic source (the Magellanic Clouds) was achieved 
by Kerr & Hindman (1953) with a 36-foot transit telescope just over half a century ago. The 
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construction of large, single dish radio telescopes produced seminal discoveries in the decade of the 
1960s, as illustrated in the fundamental paper of that period (Roberts 1975). A decade later, the 
completion of the Very Large Array (VLA) and of the Westerbork Synthesis Radio Telescope, the 
resurfacing of the Arecibo dish and rapid progress in detector and spectrometer technology made it 
possible for HI spectroscopy to achieve order of magnitude improvements in terms of sensitivity and 
resolution. New scientific problems became accessible, and extragalactic HI line research underwent 
a phase of rapid growth. The study of rotation curves led to the discovery of dark matter in spiral 
galaxies; the potential of the luminosity-lincwidth relation as a cosmological tool became apparent: 
the impact of tidal interactions and of the intracluster medium on galaxy evolution was illustrated 
in great detail through measures of the HI emission. The 21 cm line was found to be an expedient 
tool to determine accurate galaxy redshifts, playing an important role in confirming the filamentary 
nature of the large-scale structure of the Universe. The application of the luminosity-linewidth 
relation led to accurate estimates of cosmological parameters and to the characterization of the 
peculiar velocity field in the local Universe. Highly sensitive measurements in the peripheries of 
disk galaxies revealed edges in their visible components, and a number of optically inert objects 
was discovered. 

Until a few years ago, however, comprehensive wide angle surveys of the extragalactic HI sky 
were unavailable. At the close of the last decade, the advent of multifeed front-end systems at L- 
band finally made possible the efficient coverage of large sections of the extragalactic sky. The first 
such system to be used for that purpose was installed on the 64 m Parkes telescope in Australia, 
and has produced the excellent results of the HIPASS survey (Barnes et al. 2001; Meyer et al. 
2004). A second 4- feed system on the 76 m LovcU Telescope at Jodrell Bank produced the HIJASS 
(Lang et al. 2003) survey. The 1990s upgrade of the Arecibo telescope, which replaced its line 
feeds with a Gregorian subreflector system, made it possible for that telescope to host feed arrays, 
as proposed by Kildal et al. (1993). Eventually built and installed at Arecibo in 2004, this 7-beam 
radio "camera", named ALFA (Arecibo L-band Feed Array), is now operational, enabling large- 
scale mapping projects with the great sensitivity of the 305-m telescope. A diverse set of mapping 
projects are now underway, ranging from extragalactic HI line, to Galactic line and continuum, to 
pulsar searches. Here, we introduce one of these newly-initiated surveys, specifically designed to 
map approximately one fifth of the sky in the HI line, out to a distance of 250 Mpc. The survey, 
currently underway at Arecibo, is referred to as ALFALFA, the Arecibo Legacy Fast ALFA Survey. 

As illustrated in Figure 1, ALFALFA aims to cover 7074 deg^ of the high galactic latitude sky 
between 0° and 36° in Declination, requiring a total of 4130 hours of telescope time. Exploiting the 
large collecting area of the Arecibo antenna and its relatively small beam size (~ 3.5'), ALFALFA 
will be nearly eight times more sensitive than HIPASS with ~four times better angular resolution. 
Furthermore, its spectral backend provides 3 times better spectral resolution (5.3 km s~^ at z = 0) 
over 1.4 times more bandwidth. These advantages, in combination with a simple observing tech- 
nique designed to yield excellent baseline characteristics, flux calibration and HI signal verification, 
offer new opportunities to explore the extragalactic HI sky. Data taking for ALFALFA was initiated 
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in February 2005, and, in the practical context of time allocation at a widely used, multidisciplinary 
national facility like Arecibo, completion of the full survey is projected to require 5-6 years. 

As discussed in Section 4.2, simulations predict that ALFALFA will detect some 20,000 ex- 
tragalactic HI line sources, from very nearby low mass dwarfs to massive spirals at z ~ 0.06. The 
survey is designed specifically to determine robustly the faint end of the HI mass function (HIMF) 
in the local universe at masses Mhi < 10® Mq, and will at the same time provide a census of 
HI in the surveyed sky area, making it especially useful in synergy with other wide area surveys 
such as SDSS, 2MASS, GALEX, ASTRO-F, etc. In conjunction with optical studies of compara- 
ble volumes, ALFALFA will help determine the true census of low mass satellites and the widely 
distributed dwarf galaxy population in the Local and surrounding groups. Its dataset will also 
provide the basis for studies of the dynamics of galaxies within the Local and nearby superclusters, 
will allow measurement of the HI diameter function, and will enable a first wide-area blind search 
for local HI tidal features, HI absorbers at z < 0.06 and OH megamasers in the redshift range 
0.16 < z < 0.25. Survey details and status can be found by visiting its website^. 

Survey efforts of this scale and scope require careful optimization of their operational strategy 
towards achieving the science objectives within the constraints imposed by practical observing 
conditions and requirements. In this paper, we introduce the science objectives of ALFALFA, 
the principal constraints which set its strategy and the results of survey simulations which allow 
prediction of its eventual results. In a companion paper (Giovanelli et al. 2005; Paper II), we 
present results obtained during a precursor observing run, designed to allow us to test and optimize 
the ALFALFA strategy during the ALFA commissioning phase in fall 2004. 

We summarize, in Section 2, the main scientific motivations of the survey. Technical details 
of the hardware are given in Section 3, while criteria leading to the design of the survey, in the 
form of scaling laws and survey simulations are described in Section 4. Observing modes, sky tiling 
and data processing plans are presented in Section 5.3, while Section 6 summarizes sensitivity 
numbers at various stages of the survey. We elaborate on the treatment of candidate detections 
and follow-up observations in Section 7 and summarize in Section 8. Throughout the paper, we 
assume Ho = 70 km s~^ Mpc~^. 

2. Overview of ALFALFA Science Goals 

Following on the results of a number of important and sucessful previous blind HI surveys, 
extragalactic HI surveys with ALFA will exploit Arecibo's huge collecting area to explore larger 
volumes of the universe with greater sensitivity and higher angular and spectral resolution if they 
are to breaJi new science ground. With Arecibo's tremendous sensitivity and beam size advantages, 
ALFALFA is designed for wide areal coverage, thereby increasing the volume sampled locally. 



^ http:// egg. astro. Cornell, edu/alfalfa 
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yielding a deep, precise census of HI in the local Universe to the lowest HI masses. 

ALFALFA aims to survey 7074 deg^ of sky at high galactic latitudes which lie within the 
declination limits of the Arecibo telescope, —2° < Decl. < +38°, as illustrated in Figure 1. Based 
on simulations described in Section 4.2 and now verified by the results of the ALFALFA precursor 
observations presented in Paper II, ALFALFA is expected to yield on the order of 20,000 HI line 
detections, sampling a wide range of sources from local, very low HI mass dwarfs to gas-rich massive 
galaxies seen to z ~ 0.06 (~250 Mpc). HI spectra provide redshifts, HI masses and rotational widths 
for normal galaxies, trace the history of tidal events with high kinematical accuracy and provide 
quantitative measures of the potential for future star formation via comparative HI contents. As a 
blind HI survey, ALFALFA will not be biased towards the high surface brightness galaxies typically 
found in optical galaxy catalogs and moreover, in contrast to HIPASS and HIJASS, will have 
adequate angular and spectral resolution to be used on its own, generally without the need for 
follow-up observations to determine identifications, positions and, in many cases, HI sizes. The 
wide areal coverage of ALFALFA overlaps with several other major surveys, most notably the 
Sloan Digital Sky Survey (SDSS), 2MASS and the NVSS. The catalog products of ALFALFA will 
be invaluable for multiwavelength data mining for a wide spectrum of purposes, and a key clement 
of our overall collaborative program is to provide broad application, legacy data products that will 
maximize the science fallout of the ALFALFA survey. 

A primary objective of ALFALFA is the robust determination of the faint end of the HI mass 
function (HIMF). The HIMF is the cosmic number density, per bin of HI mass, of detectable HI 
line signals in a survey sensitive to the global neutral hydrogen within a system. The most recent 
estimates of the HIMF based on significant numbers of galaxies have been presented by Zwaan et 
al. (1997, hereafter Z97), Rosenberg & Schneider (2002; hereafter RS02), Zwaan et al. (2004, 
hereafter Z04), Zwaan et al. (2005, hereafter Z05) and Springob et al. (2005a). The Z04 and 
Z05 HIMFs are based on the HIPASS survey, while the RS02 and Z97 HIMFs are both based 
on surveys conducted at Arecibo during the period of its recent upgrade. The faint end slope of 
those determinations of the HIMF vary between —1.20 and —1.53, yielding extrapolations below 
Mhi = 10^ Mq that disagree by an order of magnitude near 10^ M©, the RS02 HIMF having the 
steeper slope. All of the previous HI blind surveys sample a lower mass limit just below M^/ = 10^ 
Mq. No extragalactic HI sources were detected by RS02 or Z97 with Mhi < Mq, while 3 
are claimed by Z04, and only a small number of detections have Mhi < 10^ Mq. We note that 
the distances of those detections are highly uncertain, for they are very nearby and the impact of 
peculiar velocity on the observed redshift is quite large, as pointed out by Masters et al. (2004). 
Thus current inferences on the behavior of the HIMF at low mass levels are quite unreliable, as 
they are based on very few objects of highly uncertain distance. 

With the aim of exploring the HIMF at masses Mhi < 10^ Mq, ALFALFA will cover a very 
large solid angle in order to survey an adequate volume at D < 20 Mpc, a distance within which 
the low HI mass systems are detectable. As shown by the simulations described in Section 4.2, 
ALFALFA will detect several hundred objects with Mhi < 10^'^ Mq. In addition, its extensive 
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catalog of more massive objects will allow comparison of the high mass end, Mhi > 10^ Mq, in the 
diverse range of environments found in the volume out to 250 Mpc. The Arecibo sky to be surveyed 
by ALFALFA, as shown in Figure 1 includes the rich central regions of the Local Supercluster and 
the nearby low density anti- Virgo region as well as a number of more distant large scale features, 
most notably the main ridge of the Pisces-Perseus Supercluster and the Great Wall connecting the 
Abell 1367 - Coma and Hercules superclusters. 

At the Virgo distance, ALFALFA should detect galaxies with HI masses as low as Mq ~ 10^ 
Mq. alfalfa will cover more than one thousand deg^ around Virgo, yielding a database of 
unprecedented breadth for combination with the SDSS, GALEX and other surveys to construct a 
complete census of baryon bearing objects in the cluster and its full infall region. The combination 
of HI content, HI distribution and the derived kinematical information with other multiwavelength 
studies will enable detailed modelling of the relative efficiency of gas stripping mechanisms such as 
tides, ram pressure or galaxy harassment as the origin of gas deficiency in Virgo. ALFALFA's HI 
maps will trace intriguing HI features like the Virgo "dark cloud" (Davies et al. 2004; Minchin 
et al. 2005), the HI "plume" around NGC 4388 (Oosterloo & van Gorkom 2005) and the huge 
envelope surrounding NGC 4532 and DD0137 (Hoffman et al. 1992). In more quiescent regions 
than Virgo, extensive tidal features such as the Leo Triplet (Haynes, Giovanelli &: Roberts 1979), 
and enigmatic systems such as the 200 kpc "Leo ring" (Schneider et al. 1983) may be found. 
ALFALFA will enable the first truly blind survey for HI tidal remnants with both sufficient angular 
resolution and wide areal coverage to verify their nature. 

While HI appendages uncover past disruptive events in galaxy evolution, extended gas disks 
around galaxies represent a reservoir for future star formation activity. In contrast to HIPASS 
and HIJASS which were limited by much poorer angular resolution (15.5' and 12', respectively), 
the 3.5' beam of ALFA will resolve the HI disks of ~500 gas-rich galaxies, allowing a quantitative 
measure of their HI sizes (Hewitt et al. 1984) and the derivation of the HI diameter function. 
In combination with optical photometry, ALFALFA will determine the fraction of galaxies with 
extended gas disks and enable studies of their host galaxies, their environments, morphologies and 
the role of gas in their evolution. More extremely extended gas disks, such as those found in 
DDO 154 (Krumm k Burstein 1984), UGC 5288 (van Zee 2004) and NGC 3741 (Begum et al. 
2005) may lurk yet unidentified. Because of its wide sky coverage, ALFALFA will trace important 
high- velocity cloud (HVC) structures in and around the Milky Way, such as the northern portions of 
the Magellanic Stream and Complex C at several times better spatial and spectral resolution than 
HIPASS, particularly important advantages in the case of narrow linewidth HVC cores (Giovanelli 
Sz Brown 1973). Because of its high flux sensitivity, ALFALFA will be eight times more sensitive 
than HIPASS to unresolved small clouds, or ultra-compact HVCs. While Arecibo cannot reach as 
far north as M31, ALFALFA will cover part of the region containing the clouds in its periphery 
identified by Thilker et al. (2004) and their possible extension toward the region around M33 
(Westmcier, Braun & Thilker 2005). 

In addition to the study of HI in emission, ALFALFA will provide a dataset well-suited for a 
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blind survey of HI absorption out to 2: 0.06. The background continuum source counts in the 
ALFALFA survey region at 1.4 GHz yield over 2000 sources brighter than 0.4 Jy and more than 
10000 brighter than 0.1 Jy. The major practical difficulty with HI line absorption studies is spectral 
baseline determination in the presence of standing waves. The large number of continuum sources 
present in the ALFALFA dataset and the adopted "drift" technique (Section 5.1) will aid in the 
assessment of whether a given spectral feature is real absorption. 

By its combination of studies of HI emission and absorption in the local universe, ALFALFA 
will allow a robust estimate of the local HI cross section, as well as a measure of its clustering 
correlation amplitude and scale. 

In addition to HI line studies, the frequency range of the ALFALFA survey will also include, 
serendipitously, lines from OH Megamasers (OHM) arising from the nuclear molecular regions 
in merging galaxy systems. Approximately 100 such sources are known to date, half of which 
were discovered recently at Arccibo (e.g. Darling & Giovanclli 2002). Observations of OHMs 
hold the potential for tracing the merger history of the Universe since the sources are associated 
with merging galaxies. An essential tool in this exercise is the OHM luminosity funtion at low 
z. ALFALFA should detect several additional dozen OHMs in the redshift interval 0.16-0.25, and 
allow a more robust determination of the low z OHM luminosity function than currently available. 

3. ALFA: The Arecibo L— Band Feed Array and its Spectral Line Backend 

The construction of the Gregorian subrefiector system for the Arecibo telescope, completed 
in the late 1990s, made possible the development of focal plane feed arrays (effectively, creating a 
focal plane). This development was foreseen during the planning phases of the Gregorian upgrade 
(Kildal et al. 1993). A seven feed array was commissioned at the Observatory during 2004. Six 
of the seven feeds (numbered 1 through 6) are physically arranged on the corners of a regular 
hexagon, while the seventh (feed 0) is at its center, as shown in Figure 2. The feeds can receive 
dual, linear polarizations and their spectral response is optimized for the range 1225-1525 MHz. 
They arc stepped TEn mode horns of 25 cm aperture, as described in Cortes-Medellm (2002). 
Because the optical design of the Gregorian subreflectors maximizes the illuminated area of the 
primary by sacrificing its circular symmetry (the illuminated area is elliptical), a circular pattern 
in the sky maps on the focal plane as an ellipse of axial ratio 1.15; reciprocally, the footprint of 
the centers of the outer beams of the ALFA array on the sky is that of a hexagon inscribed in an 
ellipse of that axial ratio. Similarly, the seven beams have an elliptical shape of the same axial ratio 
and orientation as the array pattern. The major axis of the ellipse is linked to the azimuth of the 
receiver, so its orientation on the sky changes with telescope configuration. In Figure 2 the relative 
location of the beams is shown when the array is positioned at the meridian and rotated about 
its symmetry axis by 19°. In this sketch, the outlines of the beams are shown at the half-power 
response, for which the beam sizes are 3.3' along the azimuth direction and 3.8' along the zenith 
angle directon, with small variations from one beam to the other. The central beam has higher 
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gain (~ 11 K/Jy) than the peripheral beams 1-6 (gain of ~ 8.5 K/Jy), which is illustrated in 
the sketch by the brighter contours. The dotted lines indicate the tracks of constant Declination 
made by each of the beams, when data is acquired in drift mode. Projected on the sky, the ALFA 
footprint in this configuration is such that beam 1 points farthest to the North and beam 2 farthest 
to the West, for observations South of the Zenith. For observations North of the Zenith, beam 1 
points farthest to the South and beam 2 farthest to the East. 

Figures 3 and 4 show the pattern for each of the ALFA beams, obtained by mapping the 
radio source 3C 138 near transit. Sidelobe levels are very different for each of the beams of ALFA. 
Located at the center of the array, beam has the most symmetric beam pattern, with a first 
sidelobe ring near 15 dB below the response at beam center, as shown in Figure 3. Contour levels 
are plotted at intervals of 3 dB. The outer beams have a very marked comatic aberration, as shown 
in Figure 4. The first sidelobe ring of the outer beams is strongly asymmetric, reaching levels near 
7-8 dB below peak response, on the section away from the array center. This feature of the system 
will require careful attention, especially in the analysis of data obtained in the vicinity of strong 
and/or extended sources. 

The system temperature ranges between 26 and 30 K for all beam/polarization channels, when 
pointing away from strong continuum sources. 

The array can be rotated about its axis, centered on beam 0, and thus the relative position 
of the beams on the sky can be rotated along the elliptical perimeter. In the case of drift mode 
observations, it is desirable to position the array in such a manner that the beam tracks are equally 
spaced in Declination. Because of the ellipticity of the array pattern on the sky, the separation 
between beam tracks depends on both the array rotation angle as well as on the array azimuth. 
When the telescope feed arm is stationed along the local meridian, the optimal array rotation angle 
is 19°, as shown in Figure 2. In that case, beam tracks are spaced 2.1' in Declination. A single 
drift will thus sweep seven nearly equidistant tracks covering 14.5' in Declination, at slightly below 
the Nyquist sampling rate. ALFALFA will map most of the extragalactic sky in drift mode with 
ALFA stationed along the local meridian, at a local azimuth of cither 0° (for observations South of 
Zenith) or 180° (for North of Zenith). Only Declination tracks transiting within 2° of the Zenith 
will be mapped with ALFA at an azimuth near 90° or 270° in order to avoid impractically small 
zenith angles. The beam separation for equidistant tracks is 1.8' at these azimuths, and thus the 
Declination sampling will be denser, as the elliptical pattern of the sky footprint of both the array 
and the individual beams will have its major axis oriented nearly parallel to the drift direction. 
Our survey will thus be done with ALFA in only two sets of configurations: one for all observations 
between Declinations 0° and 16°, as well as between 20° and 36°, with ALFA on the meridian, and 
a second for observations between Declinations 16° and 20°. 

Spectra will be recorded every second, yielding approximately 14 samples per beam in the Right 
Ascension direction. This sampling rate, which largely exceeds Nyquist, is principally motivated by 
the advantages deriving in the identification of radio frequency interference (RFI). Further details 
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on ALFA can be found at the NAIC website ^. 

ALFALFA uses the feed array connected through a fiber optics IF-LO system to a spectral 
Une, digital backend consisting of a set of processors each individually referred to as a WAPP 
(Wideband Arecibo Pulsar Processor). The full spectral backend consists of four WAPP units, 
each capable of processing the two polarization signals from two ALFA beams. The WAPP set can 
thus produce 16 autocorrelation spectra, each with a maximum bandwidth of 100 MHz, over 4096 
lags. Fourteen of those are matched to the seven polarization pairs from the ALFA beams, and 
a spare pair duplicates the signal of the seventh beam. At the offline processing stage, the extra 
pair of spectra are used for RFI monitoring purposes. Each data record thus consists of 65,536 
spectral samples (16 x 4096). Since, as mentioned above, ALFALFA records data every second, 
the generation of raw data by the survey is slightly over 1 GB per hour, including headers. 



4. Survey Design 

The strategy for the ALFALFA survey has been developed over the last few years, balancing 
the practical realities involved in using the Arecibo telescope, the constraints of telescope time 
availability, and the principal science objectives outlined in Section 2. Here we review the consid- 
erations that enter into the survey design and numerical simulations that have been used to refine 
it. 



4.1. Scaling Relations 

The HI mass of an optically thin HI source at distance Dmpc, in solar units, is 

Mhi/Mq = 2.356 X lO'Dl^p^J S{V)dV, (1) 
where S{V) is the HI line profile in Jy and V is the Doppler velocity in km . To first order, 

Mhi/Mq ~ 2.4 X 01^,3 peakWkms, (2) 

where Speak is the line peak flux and Wtms its velocity width in km s^^ . For detection, the signal- 
to-noise ratio s = fa Speak / S noise Hiust cxceed some threshold value; /a < 1 quantifies the fraction 
of the source flux detected by the telescope's beam. The parameter fp = 1 for a point source, 
while for resolved sources, it decreases roughly like the ratio between the beam solid angle and the 
solid angle subtended by the source. An estimate of Sndse can be obtained from the radiometer 
equation for the rms figure 

Q _ O^sys/G) , , 

i~>rms — /n T— f T Ti \ ) 

V2 X A/c/j X is X ft 



^ http:/ / alfa. naic. edu 
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where Tgys/G is the system temperature divided by the system gain (for the ALFA feeds, Tgys/G 
will vary between 2.65 and 3.40 Jy; here we adopt a flat value of 3.25 Jy); Af^h is the channel 
bandwidth in Hz and tg the integration time in seconds. The factor 2 under the square root indicates 
that two independent polarization channels are averaged. For ALFALFA, Afch = 25 kHz, which at 
the rest frequency of the HI line, is equivalent to 5.3 km s"^ . The factor ft accounts for spectral 
smoothing of the signal, fsmo, the switching technique applied for bandpass subtraction, f switch, 
and other observational details, such as autocorrelation clipping losses, i.e. ft = f switch fsmo father ■ 
For the data taking scheme of ALFALFA, fswitchf other — 0.7. The signal-to-noisc of a feature of 
width Wkms is best rendered when the noise is measured after smoothing the signal to a spectral 
resolution on order of Wkms/'^- In practice, however, the smoothing of L-band spectra of Wkms — 
several hundred km does not reduce the noise in proportion to W^^^^ and, moreover, Sp^ak 
is depressed by such smoothing, for spectral shapes are by no means boxlike. The fact that the 
detection criterion described above applies well to narrow lines but not so to wider ones was also 
noted by Rosenberg & Schneider (2002). We assume here that spectral smoothing will increase 
signal-to noise up to a maximum Wkms — 200, and that smoothing beyond that width will be 
ineffective in increasing s. For a conservative signal-to-noisc threshold of 6, we can then write: 

where 7 = —1/2 for Wkms < 200 and 7 = —1 for W^ms > 200. By inverting, we can obtain a 
minimum detectable HI Mass 

{^)mir^ = ^-m^'Dl^c t-;''HWkms/m~' ■ (5) 

With an integration time of 30 sec per pixel solid angle (see Section 6) or 48 sec per beam solid 
angle, ALFALFA should thus detect an HI mass of 10^ M©, Wkms = 25, at a distance of ~ 6.5 
Mpc, and a source of 10^ M© and of the same width out to ~ 20 Mpc. 

It is useful to review some of the basic scaling relations relevant to the design of a survey: 

• The minimum integration time required to detect a source of HI mass Mhi and width Wkms at 
s = 6, at the distance Dmpc with ALFA is, from eqn. 5, 




i.e. the depth of a survey increases only as tj . With equality of back-ends, the ts required 

to detect a given Mhi at a given distance decreases as the square of G, i.e. as the 4th power of the 
reflector diameter. Arecibo offers a tremendous advantage because of its huge collecting area. 

• The beam of a telescope of collecting area ^4 is fi^ oc , while the maximum distance at which 
a given HI mass can be detected is Dmax G^/^. Since G cc A, the volume sampled by one beam 
to the maximum distance Dmax is ^beam 

oc r2foL>j^(j2,/3 oc A^/'^, i.e. in a fixed time, a radio telescope 
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samples a volume that scales with the reflector diameter, yielding a significant comparative 
advantage for a large aperture like Arecibo. 

• Assuming that clouds of mass Mhi are randomly distributed in space out to the maximum 
distance at which they are detectable, DmaxiMHi), the number of clouds detected by a survey 
increases linearly with the sampled volume Vsurvey = ^survey D^^^^f 3, where ^^survey is the solid 
angle mapped by the survey. We can thus increase the number of detections either by sampling a 
larger solid angle ^survey or by increasing Dmax{^Hi)- Now, the total time required to complete 
the survey is 

t survey ^ iS^ survey /^h)ts^ (7) 

where is the telescope beam. Since DmaxiMni) oc t^^ , as shown in equation 6, we can write 

^surveyi,-^ H CX: ^survey\Dmax{,^Hiy\ 0^ ^survey^s^ ^ tsurvey^s ^ ' (^) 

and inverting: 

tsurvey Vsurvey{MHl)Dmax{MHl) OC Vsurvey{M h 1)1^ ^ . (9) 

To achieve a given surveyed volume Vsurvey{MHi), once Mhi is detectable at an astrophysically 
interesting distance, it is more advantageous to maximize than to increase the depth 
of the survey DmaxiMni)- 

The scaling relations described above provide only general guidelines in the design of a survey. 
Other considerations can and will play important roles in the survey strategy. For example, the 
growing impact of RFI on HI spectroscopy dictates increased attention to signal identification and 
corroboration, recommending a survey with more than a single pass over a given region of sky, as 
we discuss in Section 5.2. The determination of specific properties of galaxies or systems may drive 
towards deeper surveys of narrow solid angle regions, as planned for other ALFA surveys with the 
Arecibo telescope, the goals and products of which will be complementary to ALFALFA. 



4.2. Survey Simulations 

The scaling relations described above dictate that ALFALFA cover a very large solid angle. In 

practice, the survey design must weigh the desire to cover a wide area with the need for sensitivity. 
An indispensable aid in the design of a survey is a thorough examination of expectations, vis-a-vis 
variance over the survey parameter space. To this end, we have carried out an extensive set of 
survey simulations to help in the design of ALFALFA and present a sample of the results in this 
section. 

The main ingredients for our survey simulation are: (i) the survey mode and sensitivity param- 
eters, deriving from the instrument configuration; (ii) an estimate of the space density of sources 
given by an adopted HIMF; (iii) an understanding of the clustering properties and deviations from 
smooth Hubble flow in the local Universe. Sensitivity considerations were presented in Section 4.1. 
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For the HIMF, we use two recent estimates which differ strongly from each other at the low mass 
end: that of Z97 and that of RS02. The more recent HIMFs by Z04 and Springob et al (2005a) 
arc bracketed by those of Z97 and RS02. Wc use a density map of the local Universe provided 
by Branchini et al. (1999), which is a density reconstruction derived from the PSCz catalog. The 
grid we used has a spacing of 0.9375 Mpc in the inner 60 Mpc, and a spacing twice that 
value between 60 and 120 Mpc, where h = 0.7; the map is smoothed with a Gaussian filter of 
a = 3.2h~^ Mpc. For distances larger than 120 h"^ Mpc, we assume a constant density. 

Sources are seeded using the density map and, separately, each of the two HIMFs. The HI gas is 
assumed to be optically thin. Because many of the sources will be resolved by the Arecibo beam, an 
estimate of HI sizes is necessary. Assuming that the HI distribution is disk-like, inclinations to the 
line of sight and linewidths need to be assigned to each source. Wc use empirical scaling relations 
obtained from our own HI survey data (Springob et al. 2005b) and Broeils &: Rhee (1997), we 
add random inclinations, realistic scatter and broad-band spectral baseline instability. With these 
recipes, we have inspected a wide grid of survey parameters in arriving at the adopted ALFALFA 
survey strategy. 

As an example. Figure 5 displays the variation with distance of the HI Mass of expected 
detections within the region of the ALFALFA sky coverage, as described in Section 5.3, adopting 
as input the two different HIMFs: the Z97 HIMF (panel a) and the RS02 (panel b). Only expected 
detections out to 150 Mpc arc shown. The simulation corresponds to an integration time ts = 30 
seconds per map pixel solid angle. The number of detections expected with RS02 is 22,200, while 
the number expected with Z97 is 15,022, with a detection threshold of S/N = 6, as defined in 
Section 4.1. The difference in expected detections is more dramatic when the HI mass of the source 
is restricted to Mhi < 10^ Mq. In that case, we expect 1400 detections with RS02 and only 249 
with Z97. It is interesting to point out that in the course of ALFALFA precursor observations, 
reported in the companion Paper II, three objects with Mfji < 10^ Mq were detected. Albeit of 
still marginal statistical value, that rate is consistent with the high end of the expectations (RS02 
HIMF) obtained from the simulations reported here. 

Three curves are inset in the panels of Figure 5: the two solid lines are the loci of constant 
integrated HI line flux of 0.9 and 1.25 Jy km s~^ , respectively. The lowest of the two corresponds to 
the completeness limit of the survey for sources of W^^ns < 200 km s~^ (detections below that line 
correspond to sources of smaller width). The second curve corresponds to completeness limit for 
sources of the same width, for an integration time per pixel of tg = 14 seconds. Such an integration 
time applies to tlic analysis of individual drift tracks, without the corroborating support (and higher 
resulting integration) of spectra in beam tracks at neighboring Declinations. This detection limit 
would result if source extraction would be carried out, for example, right after data taking, and 
before an entire data cube (spatially two-dimensional, plus one spectral dimension) is available. 
In this case, the expected number of detections would be 13,804 for the RS02 HIMF, and 9601 for 
the Z97 case, a drop of respectively 38% and 36% from the previous set of numbers. The decrease 
in the number of detections with small HI masses would be more severe if signal extraction were 
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applied to individual tracks only, rather than to full maps: in that case, only fewer than half of the 
sources would be detected in the RS02 case, and just above half in the Z97 case. 

The topmost (dashed) curve inset in Figure 5 corresponds to a flux integral of 6.8 Jy km s~^ , 
the HIPASS completeness limit at the 6a level, for detection of sources of width < 200 km s~^ ; 
this is the HIPASS analog of the lowest of the two sohd hnes for ALFALFA. It uses a HIPASS 
limit of 13.3 mJy per map pixel S^rGcL, clS reported by Barnes et al. (2001). This provides a good 
graphical illustration of the comparison between the two surveys. 

Figure 6 shows the sky distribution of the detected sources by an ALFALFA-like survey, in 
the simulation with the Z97 HIMF. Sources at all Right Ascensions are plotted, albeit ALFALFA 
will only cover 60% of the full R.A. range. In the lower panel, only the detections with Mhi < 10^ 
Mq are plotted. Figure 7 shows the analogous graphs for the RS02 HIMF. 

Simulations results such as those presented here aid us in the estimate of the statistical efficacy 
of the survey data, most importantly in the determination of the faint end of the HIMS and the 
clustering properties of the new detections. 

5. Observing Mode 

Given the science objectives outlined in Section 2 and the scaling relations and simulations 
presented in the preceding section, the final consideration of the survey design strategy takes into 
account more telescope-related practicalities. In this section, we review those issues which have led 
us to adopt a very simple observing strategy, a two-pass drift scan mode, covering the sky with 
"tiles" extending from 0° < Dec. < +36°. 

5.1. Drift Mode 

The Arecibo telescope is an altitude-azimuth system located at a latitude near 18°. Its Gre- 
gorian dome can be steered within ~ 20° of the zenith, but the system gain and performance 
degrades rapidly at zenith angles above ~ 18°. In general, the performance characteristics includ- 
ing beamwidth, pointing accuracy, sidelobe levels, spectral baseline stability and susceptibility to 
RFI vary with both azimuth, zenith angle and feed rotation angle. Furthermore, the ALFA foot- 
print on the sky and its beams' structure vary likewise in a complicated manner. The design of 
ALFA surveys is thus strongly constrained by this variance, and with it, the corresponding degree 
of calibration complexity a particular observing program can endure. ALFALFA aims to minimize 
the impact of these factors on performance through a choice of maximum simplicity in the observing 
mode and minimal electronic intrusion at the detection level. 

As mentioned previously in Section 3, the ALFALFA survey is thus being carried out in a 
fixed azimuth drift mode. For most of the survey, the azimuth arm of the telescope is stationed 
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along the local meridian, the zenith angle of ALFA determining the declination to be mapped. A 

tiny elevation readjustment is periodically applied to maintain drift tracks at constant J2000.0, 
rather than current declination. Without such adjustment, drift tracks taken few years apart 
would noticeably diverge from one another. No firing of noise calibration diodes is done during 
normal data acquisition. Rather, data taking is interrupted very briefly every 10 minutes while a 
calibration noise diode is fired for 1 second. This interruption produces data gaps of 5 seconds — 
approximately 1/3 of a beamwidth in R.A. — in each 600 second drift data stream. The central 
frequency of the bandpass is set and never changed during an observing session, i.e. no Dopplcr 
tracking of the local oscillator frequency is applied (realignment to a common heliocentric reference 
frame is applied to each spectrum off line, using high precision time, position and Earth's motion 
stamps updated every second in data headers). With no moving telescope parts, constant gain 
and nearly constant system temperature along a drift are obtained; standing waves will change 
slowly, as driven by the sidereal rate; beam characteristics remain fixed; bandpass subtraction is 
optimized. 

The solid angle mapped by the survey is subdivided into "tiles" of 4° in Declination (see below 
and Figure 1) extending from Dec.=0° to 36°. For 8 of the 9 bands of tiles, the azimuth of the 
feed array is along the local meridian (at azimuth either 0° or 180°) while the rotation angle of 
the feed array is fixed at 19°, yielding tracks for the seven array beams that are equally spaced 
in Declination. In order to map the telescope's "zone of avoidance" near zenith, the band of tiles 
centered at Dec. =+18° will require a different orientation: with the azimuth arm nearly E-W. 
This strategy greatly simplifies the disentangling of main beam and sidelobe contributions to the 
maps: characterization of ALFA parameters thus needs to be made on a greatly reduced volume 
of telescope configuration parameter space. 

Drift mode observations, combined with the calibration scheme described above, yield maxi- 
mally efficient use of telescope time, providing high photometric quality with very small overhead. 
We expect that, bar instrumental malfunctions, telescope time usage for science data will approach 
97%. 

5.2. Two— Pass Strategy 

As discussed in Section 4.1, the volume sampled at any HI mass limit, for a survey of fixed 
total duration, varies with the integration time per point as tg ■ Once a threshold sensitivity is 

reached, it is more advantagcoTis to increase the solid angle of the survey than its depth. Because of 
the spacing of the ALFA beam tracks in drift mode discussed in Section 3, coverage of the sky in a 
one-pass drift survey is slightly worse than Nyquist. For a fixed amount of observing time, a single 
pass strategy would appear to maximize the number of detected sources. For a fixed total survey 
time, the loss of survey volume sampled by going from a one-pass to a two-pass drift survey is, 
according to Section 4.1, 19%. Several advantages of a two-pass strategy offset that loss, however: 
(1) A second pass will greatly aid the separation of cosmic emission from RFI which is unlikely to 
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afFect each pass identically; (2) The denser sky sampling will allow statistical separation of spurious 
signals from cosmic ones, thus allowing reliable detection to lower values of S/N; (3) If the two 
passes arc made when the Earth is at very different phases of its orbit, confirmation of cosmic 
nature for detection candidates can be obtained by verifying that they are separated in topocentric 
radial velocity by 30cos(Af?) km s~^ , where is the change in the angle between the line of 
sight to the detection candidate and the velocity vector of Earth on its heliocentric orbit; this 
requires that the second pass be undertaken 3 to 9 months after the first pass, modulo a year; (4) 
The variability of radio continuum sources can be measured, and radio transients can be identified, 
allowing commensality with other science teams interested in studies of those phenomena; (5) Given 
design features of the ALFA hardware, maintenance will be difficult and, as a result, ALFA may 
operate at less than 100% capacity (i.e. one or more beams may be unusable) during some fraction 
of the time. Loss of a beam in a single-pass survey would result in grievous holes in sky coverage, 
whereas a two-pass strategy would greatly attenuate the resulting damage to the survey. For all 
these reasons, the high galactic latitude Arecibo sky to be mapped by ALFALFA will be covered in 
two drift passes. The resulting effective integration time of the survey, per beam area, will be about 
48 seconds. Another way of expressing the sensitivity of the survey is in terms of the integration 
time per deg^, which will be about 14,700 seconds. 



5.3. Sky Tiling and Data Products 

As shown in Figure 1, the sky to be mapped by ALFALFA extends between 0° < Dec. < 36° 
and over two blocks of Right Ascension, respectively 07'*30"' to 16'*30"' and 22'*00"' to 03'*00"', 
although the vagaries of telescope time allocation will produce some irregularities in the survey 
solid angle boundaries. The exclusion of the low galactic latitude regions within the telescope's 
horizon is driven by (a) the realistic assessment that pulsar and other galactic ALFA surveys will 
greatly increase the pressure on low galactic latitude LSTs and (b) the expectation that part of the 
low galactic latitude, extragalactic sky will be surveyed commensally with pulsar and other galactic 
surveys. 

For bookkeeping and data release purposes, the sky mapped by ALFALFA will be subdivided 
into 378 tiles, each of 20"* in R.A. and 4° in Dec. Mapping a tile in single-pass drift mode requires 
17 drifts of ALFA, spaced ~ 14' in Dec. and each yielding 7 drift tracks; equally as many additional 
drifts are required to complete the second pass at a later time. For the second pass, beam tracks 
will be interleaved with those of the first pass, so that the final Declination sampling will be ~ l', 
better than Nyquist. In order to minimize "scalloping" of the gain over the map introduced by 
the higher gain of central beam relative to the outer ones, the second pass drifts are offset by 
7'18" relative to the first pass tracks. 

The data processing environment chosen for ALFALFA is IDL. A substantial body of spectral 
line software generated by one of us (PP) already exists at the Arecibo Observatory. Further 
development specific to ALFALFA has been grafted on this fertile base. The tile size was chosen to 
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constitute a data block that can reasonably be handled for data processing in an efficient manner 
by current desktop computers. The generation of raw data proceeds at the rate of ~ 1.2 GB/hr, 
and upon conversion from its raw FITS format to an IDL structure, a single 600 sec drift is ~ 200 
MB. Such a data block is well suited for one of the most computer intensive parts of the reduction 
pipeline, that of bandpass subtraction. The data for a full tile, after polarization averaging and 
regridding, can fit within the 2-4 GB memory of current low-cost desktops. 

The data processing path for ALFALFA data can be summarized as follows: 

• One FITS file per 600-record drift is generated by the data taking software at the Arecibo 
Observatory. By the end of each observing session, each of those is converted into an IDL 
structure, a 'drift' structure, and stored for further analysis at the Observatory and the 
observers' institutions. 

• Within weeks, all data of an observing session is noise-calibrated and a bandpass solution is 
computed. The "bandpassed" , calibrated and baselined spectral data for each beam/polarization 
configuration are obtained as output of an automated pipeline that is designed to preserve 
not only small angular scale features such as external galaxies, but also large structures such 
as HVCs and galactic HI. 

• The first detailed visual inspection of the data follows, in the course of which the observer 
flags regions of each position-velocity map for RFI and other occurrences of data corruption. 
It is anticipated that as much as about half of all sources to be detected by ALFALFA will 
be visible to the eye at this stage. A first automated signal extraction algorithm pass will 
produce a list of candidate detections. Noise diode-calibrated, bandpass-corrected, baselined 
and RFI-fiagged spectra as obtained to this stage constitute what we shall refer to as Level 
I Data Products. 

• Upon completion of the second pass through a given sky tile, data will be re-calibrated 
using the continuum sources present within the tile, regridding of the sky sampling will take 
place, after smoothing by a homogenoeus resolution kernel and conversion into data cubes 
will follow. The output of this processing stage shall be referred to as Level II Data Products. 

Because telescope scheduling is a dynamic process which responds to proposal pressure at a national 
center, the scheduling of data releases far in advance is not possible. However, the ALFALFA 
observing status is continuously updated at the survey website^. Observing plans foresee completion 
by bands of tiles and, when an accurate prediction of the completion is available, data release plans 
for that band will be posted. Data release will take place through an ALFALFA/HI node connected 
to the U.S. National Virtual Observatory. A preliminary example of web-based data presentation is 
linked to the aforementioned website as well as directly reachable ^: it allows access to the spectral 



^http:/ / egg. astro. Cornell, edu/alfalfa 
*http:/ / egg. astro. Cornell, edu/precursor 
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line data, images of optical counterparts and parameter tabulation of the precursor observations' 
results discussed in Paper II. 

As an example of a Level I data product and the potential of the ALFALFA survey, Figure 

8 shows a single pass drift across the galaxy NGC 3628, one of the Leo Triplet galaxies, and its 
tidal tail. The displayed position-velocity image consists of a constant declination drift of 600 
1-sec records at Dec. (J2000) = +13°36'45", corresponding to Dec. (B1950) = +13°53'09" for 
comparison with Figure 1 of Haynes, Giovanelli & Roberts (1979). The single drift rms noise in the 
image is 3.5 mJy. Contours are linearly spaced by 6 mJy and the lowest contour is plotted at 2 mJy 
per beam. The tidal tail is traced by the ALFALFA data as far as the earlier point-by-point map, 
but a countertail at earlier Right Ascension than NGC 3628 is also clearly visible in the new map. 
After the second pass ALFALFA data become available, the Leo Triplet region will be mapped 
with a sensitivity twice as deep as that shown in Figure 8, allowing a detailed study of the system 
over a wide area. 

6. Expected Survey Sensitivity 

Here we summarize the sensitivity of ALFALFA at several survey levels: 

• A 1-second record of a drift scan, after accumulation of both polarizations, will yield a 
spectrum of Srms — 13(res/10)~-'^/^ mJy, where res is the spectral resolution in km s"^ . 

• A single drift, position-frequency map spatially smoothed to the spatial resolution of the 
telescope beam will yield Srms — 3.5(res/10)~^/^ mJy. 

• A spatially two-dimensional map of two-pass ALFALFA data, smoothed with a kernel of 2' 
at half power, will have Srms — 2.3(res/10)~^/^ mJy per pixel. 

• The rms sensitivity per beam area, after a two-pass survey, will be Srms — 1.8(res/10)~^/^. 

• The 6cj HI column density limit will be NHi,Um = 1.6 x 10i^(VF/10)(res/10)"i/2 atoms cm , 
for a spectral line of width W km s~^ , observed with a spectral resolution of res km s~^ . 

Column density sensitivity is, in general, independent of telescope size, and thus ALFALFA 
will not reach deeper Nhi levels than previous wide angle surveys such as HIPASS. In fact, given 
the shorter integration time per beam area, ALFALFA will have lower sensitivity to Nhi than did 
HIPASS for very extended sources. It may be argued that only surveys with longer integration 
times per beam area than HIPASS can break new ground. However, this argument holds true only 
if sources are well resolved by the telescope beam. The beam area of Arecibo is nearly 20 times 
smaller than that of the Parkes telescope. If sources are unresolved by the beam, the telescope 
can only detect total flux, and the observation cannot be used for any inference on source column 
density. In fact, very few extragalactic HI sources were resolved by the Parkes beam; the smaller 
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Arecibo beam size gives ALFALFA a major advantage. To illustrate this point, Figure 9 shows 
two histograms of the angular size distribution of optically selected, catalogued galaxies: the upper 
one for galaxies known to be within the ALFALFA survey region, the lower one for galaxies in the 
whole of the southern hemisphere. The optical size used for this comparison is D25 as catalogued 
in the Third Revised Catalog of Galaxies (de Vaucouleurs et al. 1991). It is well established that, 
on average, the HI size (measured at the level near 1 Mq pc~^) of optically selected galaxies is 
about 1.6 times that blue size (Broeils & Rhee 1997), although for dwarf irregular systems that 
ratio may rise significantly (Swaters et al. 2002). Even allowing for a small number of extremely 
large Hl-to-optical size ratios, the total number of galaxies resolved by the Parkes telescope beam 
over the whole of the southern sky is on the order of a dozen or two. Only for those, the Magellanic 
Stream and High Velocity Clouds is column density sensitivity of any relevance for HIPASS. The 
ALFALFA survey should, on the other hand, resolve several hundred galaxies and High Velocity 
Clouds, and map their peripherics to a column density limit of the order of 5 x 10^^ atoms cm~^. 
Through careful analysis of the HI mapping datascts including consideration of the impact of 
sidelobe contamination in the two or three selected telescope configurations adopted for the drift 
survey, the ALFALFA survey will address the issue of whether a column density regime below 10^^ 
cm~^ is commonly found in the local Universe (Corbelli & Bandiera 2002). 

7. Candidate Detections and Verification of Cosmic Signals 

ALFALFA will produce a catalog of tens of thousands of candidate detections; on order of 
20,000 will be cosmic sources. The number of candidate detections per bin of signal-to-noise s will 
increase steeply, and the probability that the candidates are real sources decreases rapidly with 
diminishing s value. 

Internal (i.e. within the survey data set) corroboration of candidate detections will rely on (a) 
comparison of independent polarization samples and (b) comparison of spatially adjacent survey 
samples. The effectiveness of part (b) will depend on the spatial sampling density and, in the 
case of multiple drifts through the same region, on the temporal consistency of the data. These 
comparisons will help exclude many marginal candidate detections of non-cosmic origin, which we 
shall refer to as "false" candidates. 

Post-survey, corroborating observations will be desirable to confirm candidate detections just 
below the signal-to-noise threshold above which signal corroboration can be internally possible. 
This may allow significant expansion of the survey 'catch' with modest additional amounts of 
telescope time. The usual compromise is necessary in setting a s threshold: too high a threshold 
will lose many valuable potential detections; too low a threshold will require impractical amounts 
of post-survey telescope time; a haphazard criterion may corrupt the completeness of acquired 
samples. 

We expect that the bulk of follow-up observations to corroborate the cosmic nature of detection 
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candidates will be carried out at Arecibo, using a single-pixel feed, hopping from one candidate 
to the next, minimizing slew and setup time. Two modes may lead to low efficiency usage of 
telescope time: too dense a set of follow-up targets may be comprised of too large a fraction of 
false candidates, and thus produce a low yield per unit of telescope time; too sparse a set may lead 
to a large fraction of the time spent slewing and in other overhead. Careful optimization will be 
required. We consider some of these issues in this section. 



7.1. Types and Numbers of Candidate Detections 

Visual inspection or automated signal extraction algorithms identify detection candidates that 
can be assigned to three classes: (i) cosmic sources, (ii) extreme statistical fluctuations of noise and 
(iii) spurious signals due to RFI or other instrumental and data analysis causes. 

Cosmic Sources. Figure 10 shows a histogram of signal-to-noise s of expected sources, as obtained 
in one of the simulations described in Section 4.2. We emphasize that the number of sources rises 
steeply as s decreases, illustrating the well-known fact that most of the survey candidate detections 
will occur near the detection limit. 

Because of the massive bulk of ALFALFA data sets, signal extraction will largely rely on 
automated procedures. Amelie Saintonge has coded a matched-filter, cross-correlation signal ex- 
traction algorithm, described elsewhere (Saintonge at al. , in preparation). Detection probability 
simulations have been carried out with this algorithm, by randomly injecting a Gaussian signal in 
simulated spectra with Gaussian noise, and monitoring the effectiveness of the signal extraction 
algorithm in recovering the injected signal. The detection probability is computed as the fraction 
of all trials in which the signal extraction algorithm positively identifies the injected signal; such 
probability is monitored as a function of s and of signal width. When the noise is measured after 
spectrally averaging over 1/2 the spectral width of the injected signal, the detection probability is 
largely independent of the signal width. 

Statistical noise fluctuations. For Gaussian noise, the probability that a single spectral channel 
yield a fluctuation of signal-to-noise between si and si + dsi is 

pidsi = ^e-^i/^dsi, (10) 

where si = Sp^ak/^^i, with Sp^ak the peak flux density and cJi the rms noise, with single-channel 
spectral resolution. Similarly, the probability for a channels-wide spectral feature to exhibit a 
deviation between s„ and s„ + dsn is 

1 2 

PndSn = -j=e~^^''^dSn, (11) 

V 27r 

where s„ = si/yn„^. 
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In a survey of Ni^g line of sight samples, taken with a spectrometer of Nc channels, the number 
of samples riw channels wide, with signal-to-noise between s and s + ds is 

ns,n^ds = Nios—Pnds = Nios—-^e~^'^/'^ds (12) 

n-w y 27r 

and the total number of statistical fluctuations of that width with s larger than a threshold Sth is 

AT 1 rCO AT 

Ns,„i = Nios — ^ / e-^ '''ds = Nios^[F{oo) - F{s)], (13) 

n-w V ivr J 

where, again, the noise is computed with a spectral resolution of n^/2 channels, and F(s) is 
the familiar cumulative distribution of the normal error function: F{—oo) = 0, F(oo) = 1 and 
F{0) = 0.5. The total number of purely statistical noise fluctuations between and sj,, with 
widths between Uwi and nw2, appearing in the survey will then be 

Nla,b],[l,2] = E Nios^iFiSa) " F{sk)]. (14) 

To first order, and ignoring the fact that in the expression above a high s, broad feature gets 
overcounted as several, lower s, narrower ones, we can approximate 

N[aM\i,2] ^ NiosN, In^ [F{sa) - F{s,)]. (15) 

For example, for a survey that samples 10^ lines of sight, with a spectrometer usefully covering 
85 MHz with Nc = 3600 spectral channels (so that a velocity width range between 25 and 500 
km s~^ translates into n«,i ~ 5 and — 100, one should expect A^>3 ~ 2 x 10^ features with 
s > 3, A'^>4 = 5 X 10^ features with s > 4, and A^>5 = 5 x 10^ features with s > 5, with width 
anywhere between 25 and 500 km s~^ . These are ominously large numbers when compared to 
expected numbers of cosmic sources between 10^ and 3 x 10^, for that range of s. 

RFI or other Spurious Signals. The above assumption of Gaussian noise is heuristic. The true 
nature of the noise will be ascertained after a significant fraction of the survey data will have 
been collected and the "normal" characteristics of both equipment and RFI environment will have 
been measured. As we discuss in Paper II, the precursor observations carried out in 2004 were 
taken in commissioning mode for the ALFA hardware, and several internal bugs have been found 
and fixed since those observations were completed. Those observations are not well suited for a 
careful analysis of the problem. For the moment, we will ignore the impact of RFI and other non- 
Gaussian sources of noise, and restrict our analysis to the discrimination between cosmic sources 
and Gaussian noise fluctuations. 

7.2. Discriminating Among Candidate Detections 

Assuming that the majority of weak detection candidates will be unresolved by the telescope 
beam, the most important means of discriminating between cosmic sources and noise fluctuations 



- 22 - 



will result from comparison of contiguous drift tracks and different polarizations of the same beam. 
Consider a single-pass ALFALFA survey, whereby contiguous drift tracks are separated by 2.1' in 
Declination on the sky. A point source swept by one of the feeds will appear also in contiguous 
tracks, at lower s. For an ALFA beam averaging 3.5' width at half power, the response 2.1' off the 
beam center is about 0.37 of that on beam center. In a two-pass drift survey, the track separations 
are 1.05', and at that distance from the beam center, the beam response is 0.78 of that on beam 
center. A point source will thus be far more easily confirmed in a two-pass survey. Similarly, since 
the HI line is unpolarized, the comparison of the two independent polarization spectra of the same 
beam will deliver equal signals, plus noise, for cosmic sources, and completely uncorrelated results 
for noise. We shall refer to the exclusion of detection candidates made possible by comparison of 
adjacent drift tracks and polarization channels as vicinity trimming, and distinguish three sets of 
detection candidates: (1) that obtained without any vicinity trimming, and those (smaller ones) 
obtained after (2) vicinity trimming in a single-pass survey and (3) vicinity trimming in a double 
pass survey. 

Figure 11 displays the cumulative number of candidate detections plotted as a function of 
signal-to-noise, expected for an ALFALFA-like survey. The simulation assumes a RS02 HIMF 
and a sensitivity corresponding to a double-pass drift survey. The thick line corresponds to the 
detection candidates associated with cosmic ('real') sources. The three thin lines correspond to 
the expected number of noise fluctuations with the given s in the three 'vicinity trimming' cases 
described above. At s > 6, most of the detection candidates are cosmic sources. At s ~ 5, the 
number of 'real' cosmic sources is 25% higher than at s = 6, but the candidate detections resulting 
from noise fluctuations is several times higher than that of real sources. Vicinity trimming can 
however drastically reduce the number of candidate detections deserving attention. Only at signal- 
to-noise levels s < 4 does the number of noise fluctuations overwhelm that of cosmic sources, after 
two-pass vicinity trimming. At this level, however, the impact of low power level RFI will play an 
important and yet quantitatively unknown role. 



7.3. Follow-up, Corroborating Observations 

Candidate detections with a comfortably high signal-to-noise threshold Sth will not need cor- 
roborating follow-up observations in order to confirm their reality as cosmic sources. Without 
considering the impact of RFI, sth may be in the vicinity of 6; consideration of the impact of RFI 
may raise Sth to higher values, in a variable manner depending on the frequency of the candidate 
signal. Candidates with s ~ sth or slightly below that level will be reobserved with the Arecibo 
telescope. To how low a level of s should re-observations be considered? A simple criterion would 
be that corroborating observations should be requested only for candidates of s such that the ex- 
pectation of confirmation, expressed in terms of detections per unit of telescope time, is at least as 
high as for the full survey. 

If a corroborating observation is to require an increase in the s from a value of, say, 5 for the 
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survey data, to about 7 for the corroborating observation, integration times of at least 1 minute per 

candidate will be necessary for corroborating observations. Observing runs to corroborate several 
hundred candidate detections at a time will thus be the norm. If the set of candidate detections to 
be checked is very sparse — say one candidate every several square degrees — slew times will be 
very substantial and bandpass-correcting observations will be required for each candidate source, 
more than doubling the required telescope time. In that case, the on-sour ce tint of order of 1 
minute may be a small fraction of the overall time required to observe each source. The Arecibo 
telescope slew times are respectively 0.4° in azimuth and 0.04° s^^ in elevation. A 1° change 
in elevation will require 25 seconds. It will thus be observationally advantageous if the sky density 
of tentative sources to be corroborated is high, e.g. on order of one per square degree or higher. 
Not only will that reduce the overhead of slew motions and settle time, but it will also allow 
for a running mean bandpass to be accumulated over a few contiguous staring observations, as 
the telescope configuration would change little between adjacent source candidates. In that case, 
allowing for slew and settle time, a corroborating observation of a single source will require on 
order of one to two minutes of telescope time. The steeply rising fraction of 'false' sources with 
decreasing s suggests that corroborating observations requesting single pixel telescope time at the 
level of approximately 10% of the request proposed for the ALFA observations will deliver optimal 
returns. 



8. Summary 

ALFALFA uses the new 7-beam Arecibo L-band feed array to carry out a wide area survey 
of the high galactic latitude sky visible from Arecibo. In addition to the all-important sensitivity 
advantage that accrues from using Arecibo, the world's most sensitive radio telescope at L-band, 
ALFA offers important and significant improvements in angular and spectral resolution over the 
available major wide area extragalactic HI line surveys such as HIPASS and HIJASS. ALFALFA is 
intended to produce an extensive database of HI spectra that will be of use to a broad community 
of investigators, including many interested in the correlative mining of multiwavelength datasets. 
It is specifically designed to probe the faint end of the HIMF in the very local Universe. 

As a result of practical considerations and simulations of survey efficiency, ALFALFA exploits 
a simple fixed-azimuth drift scanning — minimum intrusion — technique. A two-pass strategy 
will greatly aid in the rejection of spurious signals and RFI, thus minimizing the need for follow-up 
confirmation observations, evening out the scalloping in the maps that arises from unequal pixel 
gain, and offering the opportunity to use the same dataset for the statistical characterization of 
continuum transients. Initial tests of the hardware, software and survey observing mode, conducted 
in Fall 2004 during the ALFA commissioning phase as described in Paper II confirm the efficacy 
of the planned approach. The basic parameters of the ALFALFA survey can be summarized 
accordingly as follows: 
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• Sky coverage of 7074 deg^, between 0° and +36° in Declination, 7.5 to 16.5 and 22.0 to 3.0 
hrs in Right Ascension, with 3.5' spatial resolution. 

• Frequency coverage between 1335 and 1435 MHz, yielding coverage of extragalactic HI in 
redshift out to cz < 18, 000 km , with 5.3 km s~^ maximum spectral resolution. 

• Sensitivity of 1.8 x (res/10)^^/^ mJy per beam area, where res is the spectral resolution in 
km s"^ . 

• On order of 20,000 HI sources are expected to be detected by the survey. Extragalactic HI 
sources with Mhi — 10^ Mq will be detectable to a distance of 6.5 Mpc, while HI masses 
Mhi — 10^ Mq, will be detectable throughout most of the Local Supercluster, including the 
Virgo cluster and out to 20 Mpc. Several hundreds will have Mhi < 10^'^ Mq, thus allowing 
a robust determination of the faint end of the HIMF. 

• Public access data products will be produced on a continuing basis as subsets (tiles) of the 
overall survey are completed. 

Observations for ALFALFA started in February 2005, and completion of the survey is expected 
to take five to six years. Cataloguing a complete census of HI locally, pinning down the HIMF to 
the lowest masses and conducting the blind HI absorption and OH megamaser surveys will require 
completion of the full 5-year program, but even the initial 2005 allocation promises early science 
results in several important areas including: the mapping of nearly 1600 dcg^, more than 3 times 
the coverage with twice the sensitivity of the Arecibo Dual Beam Survey (Rosenberg & Schneider 
2000); a first blind census across the Virgo cluster region with a detection limit of Mhi > lO''' 
Mq at the cluster distance (assuming a width W = 30 km s"^ ); a complete search for HVCs 
around M33; the identification of gas-rich galaxies in the NGC 784 and Leo I groups; the mapping 
of the environments of 12 gas-rich galaxies with Dugc > 7'; and a first attempt at a large blind 
survey for HI absorbers. 
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Fig. 1. — Proposed sky coverage of the ALFALFA survey, in the Virgo (upper) and anti- Virgo 
(lower) directions. In each panel, the thicker lines at constant RA or Dec. outline the proposed 
survey area. Dashed lines at constant Dec. make the designated ALFALFA 'tile' strip boundaries. 
The thick dotted curves to the right of the upper panel and top of the bottom panel mark b 

= +20° (upper) and —20° (lower) while the set of three thick lines crossing each panel top to 
bottom trace SGL = —10°, 0° and +10°. Filed circles mark galaxies with observed heliocentric 
recessional velocities cz < 700 km s~^ , while open circles denote objects believed to lie within 10 
Mpc (Karachentsev et al. 2004), based largely on primary distances. 
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Fig. 2. — Sketch of the geometry of the ALFA footprint, with the array located along the local 
meridian and rotated by an angle of 19° about its axis. The outer boundary of each beam corre- 
sponds to the -3 dB level. The dashed horizontal lines represent the tracks at constant Declination 
of the seven ALFA beams, as data is acquired in drift mode. 
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Fig. 3. — Beam pattern of beam 0. Contour lines and shading intervals are plotted at intervals of 
3 dB below peak response (the highest contour is at half the peak power). The first sidelobe ring, 
with a diameter near 12', is at approximately -15 dB. 
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Fig. 4. — Beam patterns of the six peripheral ALFA beams. Contour hnes and shading intervals 
are plotted at intervals of 3 dB below peak response (the highest contour is at half the peak power) . 
Note that the sidelobe levels are significantly larger than for the central beam 0, and that they rise 
steeply on the outer side of the array, exhibiting strong comatic aberration. 
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Fig. 5. — HI Mass vs. distance plot of expected detections within D < 150 Mpc, assuming a Z97 
HIMF (panel a) and a RS02 HIMF (panel b). Calculations were made for the ALFALFA sky region 
specified in Section 5.3, tg = 30 seconds per map pixel area and a detection threshold S/N = 6. 
The solid lines in both panels indicate detection limits of 0.9 Jy km s^^ and 1.25 Jy km . 
The first would be near the completeness limit of the survey for sources of width < 200 km s~^ , 
the second near the completeness limit for objects of the same width, but with tg = 14 seconds 
integration. The dashed line corresponds to a flux density of 6.8 Jy km , a 6-sigma HIPASS 
detection limit for a 200 km s~^ wide source. 
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Fig. 6. — Sky distribution of expected detections for the ALFALFA survey, as rendered with a 
Z97 HIMF. In panel (a) sources of all HI masses are plotted, while in panel (b) only those with 
Mhi < 10^ Mq are shown. Note that the ALFALFA survey will be restricted to Right Ascensions 
07'* to 16.5^ and 22^ to 3^*, although the fuU range of R.A.s is shown in the figure. 
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Fig. 7. — Analogous display to that in Figure 6, except that the RS02 HIMF was used for the 
simulation. 



-34- 




Fig. 8. — Position-velocity map cutting across NGC3628 and its tidal tail at constant Dec(J2000) 
= +13° 36' 45". The lowest contour is at 2 mJy per beam, and other contours are spaced by 6 mJy 
per beam. See text in Section 5.3 for further details. 
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Fig. 9. — Histograms of the optical major blue diameter, D25, of galaxies larger than 1', in the 
ALFALFA survey region (upper) and the whole southern hemisphere. Aperture synthesis studies 
have shown that the diameter of the HI disk for optically selected galaxies is on average 1.6 times 
larger than the optical size. The Parkes 15' beam thus resolves on order of 15 galaxies over the 
whole southern hemisphere; the ALFA ~3.5' beam should resolve several hundred galaxies over the 
ALFALFA survey region. The bin size is 0.1' in both histograms. 
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Fig. 10.— S/N histogram of a 12s ALFALFA survey using the RS02 HIM. S/N bins have width 
0.1 in S/N. S/N is defined as the peak signal flux to the rms, computed in matched filter mode, 
over a spectral resolution equal to 1/2 the signal width. Only tentative detections with S/N> 2.5 
are plotted. 
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Fig. 11. — Cumulative number of candidate detections as a function of S/N, expected for an all- 
Arecibo sky drift survey. The lower (nearly flat) line refers to the 'real' sources (assuming an RS02 
HIMF) while the upper curves refer to the three cases of 'vicinity trimming' described in Section 
7.2. 



